Shuttling heat across ID homogenous nonhnear lattices with a Brownian heat motor 
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We investigate directed thermal heat flux across ID homogenous nonlinear lattices when no net 
thermal bias is present on average. A nonlinear lattice of Fermi-Pasta-Ulam-type or Lennard- Jones- 
type system is connected at both ends to thermal baths which are held at the same temperature 
on temporal average. We study two different modulations of the heat bath temperatures, namely: 
(i) a symmetric, harmonic ac-driving of temperature of one heat bath only and (ii) a harmonic 
mixing drive of temperature acting on both heat baths. While for case (i) an adiabatic result for 
the net heat transport can be derived in terms of the temperature dependent heat conductivity 
of the nonlinear lattice a similar such transport approach fails for the harmonic mixing case (ii). 
Then, for case (ii), not even the sign of the resulting Brownian motion induced heat flux can be 
predicted a priori. A non-vanishing heat flux (including a non-adiabatic reversal of flux) is detected 
which is the result of an induced dynamical symmetry breaking mechanism in conjunction with the 
nonlinearity of the lattice dynamics. Computer simulations demonstrate that the heat flux is robust 
against an increase of lattice sizes. The observed ratchet effect for such directed heat currents is 
quite sizable for our studied class of homogenous nonlinear lattice structures, thereby making this 
setup accessible for experimental implementation and verification. 

PACS numbers: 05.40-a,07.20.Pc,05.90.-|-m,44.90.+c,85.90-|-h 



I. INTRODUCTION 

In recent years, we have experienced a wealth of theo- 
retical and experimental activities in the field of phonon- 
ics, the science and engineering of phonons Ij. Tradi- 
tionally being regarded as a nuisance, phonons are found 
to be able to carry and process information as well as 
electrons do. The control and manipulation of phonons 
manifest itself in the form of theoretically designed ther- 
mal device models such as thermal diodes [1, 
thermal transistors Q, thermal logic gates 0] and ther- 
mal memories [l^. The theoretical research has been 
accompanied by pioneering experimental efforts. In par- 
ticular, the first realization of solid-state thermal diode 
has been put forward with help of asymmetric nanotubes 
prj . Owing to the transport of phonons, the heat flow 
can be controlled the same way as electric currents. 

Dwelling on ideas from the field of Brownian motors 
[13, [Hj 0, [Hj Ell, ~ originally devised for Brownian 
particle transport, - a thermal ratchet based on a non- 
linear lattice setup has been proposed in Ref. In 
absence of any stationary non-equilibrium bias, a non- 
vanishing net heat flow can be induced by non-biased, 
temporally alternating bath temperatures combined with 
a non-homogenous coupled nonlinear lattice structure. In 
the similar smrit of pumping heat on the molecular scale 
[ll,IIi,[20,|2l|,[ll, the heat flow can be directed against 
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FIG. 1: (Color-online) Schematic setup of a ID homogenous 
nonlinear lattice, being coupled to two heat baths with peri- 
odically varying temperatures at TL{t) and Tnit). 



an external thermal bias. 

With this work we propose a superior, easy to imple- 
ment Brownian heat motor that can induce finite net 
heat flow for a homogenous^ intrinsically symmetric lat- 
tice structure as depicted in Fig. [TJ In doing so, the lat- 
tice system is brought into contact with two heat baths 
which both may be subjected to time-periodic tempera- 
tures, see Fig. [TJ To obtain a finite directed heat current 
then requires a symmetry breaking. In this work we shall 
investigate two mechanisms of symmetry breaking by use 
of a temporal modulation of temperatures of heat baths. 

Section II introduces our model for directing heat cur- 
rent through ID lattice chains with temperature mod- 
ulations applied to connecting heat baths. Section III 
presents analytical adiabatic theory and extensive nu- 
merical results for the case that temperature is time- 
modulated in one bath only, i.e. our case (i). In this 
case, an intrinsic nonlinear temperature dependence of 
the heat conductivity is sufficient to induce a shuttling 
of heat. Our main results are presented with Section IV: 
A more intriguing mechanism comes into play when ap- 
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plying an unbiased modulation of temperatures in both 
heat baths. Put differently, in order to eliminate the pos- 
sibility of solely creating a non-vanishing net heat flux 
from the the temperature-dependent thermal conductiv- 
ity k(T) wc invoke a nonlinear harmonic mixing drive of 
temperature in both heat baths; i.e. our case (ii). The 
results are discussed and summarized in the Conclusions. 



II. SHUTTLING HEAT DESPITE VANISHING 
AVERAGE TEMPERATURE BIAS 

Explicitly, we study numerically a ID homogenous 
nonlinear lattice consisting of N atoms of identical mass 
m. This setup is depicted in Fig. [U with the following 
nonlinear lattice Hamiltonian: 

N 2 ^ 
i=l 1=0 

where Xi is the coordinate for i-th atom, and the 
distance between two neighboring atoms at equilibrium 
gives the lattice constant a. The interaction between 
nearest neighbors is described by V{xi+i,Xi). Here, fixed 
boundary conditions xq = and a;Ar+i = (N + l)a have 
been employed. The first and last atom are put into con- 
tact with two Langevin heat baths, generally possessing 
time-dependent temperatures TL{t) and Tii{t), respec- 
tively. Moreover, Gaussian thermal white noises obeying 
the fluctuation-dissipation relation are used; i.e., 

(ei(w)W) = 

and (ei(JV)(i)ei(Ar)(0)) = 2kBvTLiB.)S{t) . (2) 

Here, fc^ denotes the Boltzmann constant and rj is the 
coupling strength between the system and the heat bath. 

The time- varying heat bath temperature Ti(<) and 
Tii{t) are chosen as periodic functions Ti^i^fi^{t) = 

TL(R){i + 27r/w) where To = 7Y(t) = Tfl(i) is the tem- 
porally averaged, constant environmental reference tem- 
perature. Clearly, the coherent changes of these bath 
temperatures occur on a time-scale that is much smaller 
than the (white) thermal fluctuations itself. Importantly 
in present context, this so driven system dynamics ex- 
hibits a vanishing average thermal bias; i.e.. 



/^T[t)=TL{t)~TR[t) = Q . (3) 

The time-dependent; asymptotic heat flux Ji{t) is as- 
suming the periodicity of the external driving period 
27: /lj after the transient behavior has died out This 
feature is conflrmcd in our numerical simulations (not 
shown here) . At those asymptotic long times, the result- 
in^ heat flux equals the thermal Brownian noise average 

Ji{t) = {x^^V{x,+l,Xi)/^Xi) . (4) 



The stationary heat flux J then follows as the cycle av- 
erage over a full temporal period; i.e., 



J(t) := J ^ — {xidV{xi+i,Xi)/dxi) dt 

27r Jo 

uj /■^''/'^ 

= 2^]^ Mt)dt, (5) 

which after averaging becomes independent of atom posi- 
tion i. With ergodicity being obeyed, this double-average 
equals as well the long time average 

J ^ XidV{xt+i,Xt)/dxt 
1 /■* 

= hm - / Xt{t)dV{xi+i,Xi)/dxi\tdt . (6) 

Here, the temporal average is over corresponding stochas- 
tic trajectories entering the relation in Eq. We em- 
phasize that the resulting ratcheting heat flux involves 
an average over the Brownian thermal noise forces. Put 
differently, the flux is not determined by a determinis- 
tic molecular dynamics but rather by the driven nonlin- 
ear Langevin dynamics following from the nonlinear lat- 
tice dynamics in Eq. ((T|) and being complemented with 
Langevin forces obeying the fluctuation-dissipation rela- 
tion in Eq. This in turn defines our Brownian heat 
motor dynamics. 

III. ROCKING TEMPERATURE OF ONE 
HEAT BATH ONLY 

Wc start out by considering that only the tempera- 
ture of the left heat bath is subjected to a time-varying 
modulation, i.e. our case (i) is defined by setting 

TL{t):=TL = Toil + Ai cos (ujt)] 

Tnit) Tr = To . (7) 

The symmetric temperature difference AT{t)/2To over 
one period is depicted in Fig. HJa) as the dashed line. 

The Fermi-Pasta- Ulam /3 (FPU-/?) lattice is used for il- 
lustration, where the potential term ^(3:^+1,2;^) assumes 
the following form: 

k (3 
V{xi+i,Xi) = -{xi+i-Xi-af + -{x^+i-Xi-aY (8) 

We introduce the dimensionless parameters by measuring 
positions in units of a, momenta in units of [a(TOfc)^''^], 
temperature in units of [ka?/kB\, spring constants in 
units of fc, frequencies in units of [(fc/m)^/^] and energies 
in units of [ka?]. The equations of motion are integrated 
by the symplectic velocity Vcrlet algorithm with a small 
time step h = 0.005. The system is simulated at least 
for a total time Uot = 2 • 10®. The chosen optimal cou- 
pling strength of the heat bath is fixed at 77 = 0.5. In 
the following simulations, we will take the dimensionless 
parameters k = 1, (3 = \. 
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FIG. 2: (Color-online) Modulating the bath temperature. 
The temperature bias AT(t) = Thit) — Tnit) is depicted over 
a full driving period without and with a second harmonic 
driving term, panel(a). The strength Ai = 0.6 is chosen for 
both cases. The third moment of the harmonic mixing signal 
{AT{t)/2To)^ is shown with panel (b). Note that in distinct 
contrast to the unbiased first moment the cycle average of this 
odd third moment is now nonvanishing. 



The net heat flux J as a function of driving frequency uj 
is depicted in Fig. [3l In the adiabatic hmit a; ^ 0, a neg- 
ative ratchet heat flux is obtained. In the high frequency 
hmit w — > oo, the left- and right-end atoms wiU experi- 
ence a time-averaged constant temperature. This corre- 
sponds to effective thermal equilibrium, yielding J — > 
when to —^ oo. In the non- adiabatic limit we even observe 
a reversal of heat flux, although not very pronounced. 

When the temperature is modulated slowly enough, 
the periodic influence of dynamical thermal bias can 
be viewed as the average, integrated quasi-stationary 
flux induced by the momentarily static thermal bias, i.e. 
Ji{t) = Ji{AT{t)/N). In the linear transport regime, this 
flux can be expressed as a linear transport Law, reading 
J.^{t) = K {{TL{t) + Tr))/2) AT{t)/N. Here the tempera- 
ture value of the size (N)-dependent thermal conductivity 
of the FPU lattice must be determined numerically at the 
midpoint of the, in this adiabatic case, linearly varying 
temperature profile Thus, the adiabatic net heat 

fiux assumes in leading order the result 
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FIG. 3: (Color-online) One-sided temperature modulation. 
The cycle averaged ratcheting heat flux J in a nonlinear FPU- 
/3 lattice is depicted as the function of the driving frequency 
uj. The hollow (red) circles are the results for harmonic drive 
Ai cos{ujt) with Ai = 0.5; the filled (blue) circles are the 
results for a weaker driving strength at Ai = 0.2. The simu- 
lations are performed on a FPU-/3 lattice with A'' = 50 atoms 
and with a reference temperature set at To = 0.5. The arrows 
pointing towards the left mark the adiabatic linear transport 
result from Eq. see text. Inset is the figure for Jad as the 
function of Ai and Jad has the basic units [10~^1. 
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(9) 



The predicted sign on last line in the above ex- 
pression originates from the following reasoning: For 
the considered FPU-/? lattice the thermal conductiv- 
ity possesses a temperature-dependent behavior k{T) cx 
1/T in the regime of dimensionless temperature 
T{t) < 1 [23|. With the reference temperature 
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Ai sin (ujt) 



Ai sin {tot) 



Tr 



Tr 



is alwaysf less than 

in the time window of [0,7r/Li;]. 

Therefore, a negative ratchet heat flux will result from 
the temperature modulation of only one heat bath as in 
Eq. ([7]). This prediction is corroborated with our nu- 
merical calculations in Fig. [3] for low driving frequencies 
UJ. Furthermore, the adiabatic value of Jad can be ap- 
proximately calculated from the quadrature in Eq. 
This calculated adiabatic value is marked with the ar- 
row pointing towards the left axis in Fig. [31 For small 
rocking strength, Ai ^ 1, this adiabatic prediction re- 
markably well agrees with the full nonlinear, asymptotic 
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value, see in Fig. [3l For larger driving strengths Ai the 
numerically precise full adiabatic result exhibits notable 
deviations from this linear transport Law estimate. 

Taking into account the size dependent property of 
heat conductivity k{T,N), we can express Jad further 
by taking a Taylor expansion for k (^Tq + jj^^ ]\j 

Ai sin (ciJt) 
2 



and K [Tq 
To: 



To,7V^ at reference temperature 



Jad = TT i 
ZTT 



dK{T,N) 



dT 



AiTosin {ujt) + 0{Al) 



xAiTosin {ujt)/N 
4^i^o dT ^ 

J-o 

oc A\N°'-^ 



(10) 



where in the 

dK{T,N) 
dT 



last line we have used the fact that 
assumes the same size-dependence as k(T, iV), 

i.e. n{T,N) oc A^" with a < 1 The non-zero adia- 

batic ratchet heat flux is a finite size effect since it van- 
ishes in the limit N ^ oo proportional to N°'~^. 

We numerically determined the adiabatic net heat flux 
by averaging heat flux from four lowest frequencies for 
each amplitude Ai. The amplitude effect Jad oc A^ is 
verified by our numerical calculations as can be seen in 
the inset of Fig. [3] 



IV. ROCKING TEMPERATURE IN BOTH 
HEAT BATHS 

We next consider case (ii) with unbiased tempera- 
ture modulations applied to both heat baths. The time- 
varying heat bath temperature Ti(t) and Tji{t) are cho- 
sen as: 

TL{t) := Tl = To[1 + Ai cos {cot) + A2 cos {2cut + ip)] 
TR{t) := Tb = To[l - Ai cos {ujt) - A2 cos (2wt + (p)!!!) 



where Tq — T^it) = Tji{t). The overall averaged net 
temperature difference, AT(<) = TL{t) - Tfl(f) = 0, is 
again unbiased as before. The driving amplitude Ai and 
A2 are taken as positive values and the relation A1+A2 < 
1 must be fulfilled to avoid negative temperatures. Note, 
however, that in distinct contrast to the situation with 
case (i) the time-dependent average temperature; i.e.. 



Ta^{t) = [TL{t) + TR{t)]/2^To 



(12) 



is now time-independent. This feature excludes us from 
estimating the adiabatic heat flux within a linear trans- 
port Law of heat as exercised under case (i). Put dif- 
ferently, the temperature dependent, nonlinear thermal 
conductivity k(T) alone is not sufficient to set up an 
adiabatic heat flow, see below. A resulting finite heat 



flux is thus beyond the mere role of a nonlinear lattice 
and instead is the outcome of the nonlinear interplay of 
harmonic mixing of the two frequencies in the nonlinear 
lattice dynamics to yield a zero-frequency response for 
the time-averaged nonlinear heat flow as defined by Eqs. 

In more detail: The second harmonic driving 
A2 cos(2wt -|- (p) causes nonlinear frequency mixing. The 
temperature signal AT{t) notably is unbiased with zero 
average. It causes, however, a dynamical symmetry 
"2^ 27 h thus giving rise to directed transport 
' '29I" 301 ■ The time evolution of the harmonic 



breakin; 
[2i,[23,l2 

mixing signal is depicted in Fig. [2Ua) for a phase shift of 
(p = with the solid line. The second harmonic drive may 
typically contain a non-zero phase shift (p. In accordance 
with previous studies in single particle Brownian motors 
[28l . [30 1 the resulting current is expected to become max- 
imal for ip = 0. In the following we therefore stick in 
our numerical studies, if not stated otherwise, to a fixed 
vanishing phase-shift tp = 0. Note that this drive with 
</? = is symmetric under time-reversal t — > —t; never- 
theless, time reversal is broken by the frictional Langevin 
dynamics acting upon the end atoms, see above. 

The fact that a finite heat flux results when driven by 
harmonic mixing can be reasoned physically by noting 
that in contrast to the case with A2 = 0: W hen A2 ^ 0, 
despite the vanishing temporal cycle average AT{t) = 0, 
one deals with temporal temperature differences AT{t) 
that are no longer symmetric around AT{t) = 0, cf. 
[21(a). With this harmonic mixing modulation, all odd- 
numbered moments (AT(t))(2»+i) 7^ , ?i > 1 are non- 
vanishing after the temporal cycle average. With heat 
flow in nonlinear systems typically being a function of the 
temperature bias beyond linear response regime, we thus 
nevertheless expect a net finite heat flux. This feature 
follows by observing that a leading nonlinear response 
due to the no n-vanishing tem poral cycle average of the 
third moment {AT{t)/2To)^ = (3/4)^4? A2 cos tp, see in 
Fig. [H^b), is non- vanishing. In clear contrast to a single 
particle case, see in Ref. [2^, |3l|, the amplitude A2 in 
our case (ii) assumes within lL(i?,) both signs. Therefore, 
one principally cannot even predict a priori the sign of 
the resulting heat current. This very fact is the bench- 
mark of a truly Brownian heat motor where the external 
temperature modulation is only weakly coupled to the 
Brownian motion induced heat flow [ij, [l^, [ij, [3i [3 ■ 



A. Shuttling heat across a Fermi-Pasta-Ulam /3 
lattice 

We start the study of case (ii) by considering the FPU- 
f3 lattice of Eq. ([5]) by using the same parameters as 
before. The net heat flux J as a function of driving fre- 
quency uj is depicted in Fig. ID For A2 = 0, there in- 
deed emerges no finite heat flux; this corroborates with 
theory because of the absence of dynamical symmetry 
breaking between positive and negative temperature dif- 
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FIG. 4: (Color-online) Ratcheting heat flux in the FPU-/3- 
lattice. Cycle averaged heat flux J as the function of the 
driving frequency ui. The solid circles are the results for sym- 
metric harmonic driving with A2 = Q. The hollow circles are 
the results for harmonic mixing with A2 = 0.3 and a relative 
phase shift ip = Q. The harmonic driving parameter is set 
for both cases at Ai — 0.6. The simulations are performed 
on a FPU-/3 lattice with A*' = 50 atoms and with a reference 
temperature To — 0.5. The three arrows indicate the distinct 
driving frequencies used in Fig. [5] 



ferences AT{t), see Fig. (3] The non-zero second har- 
monic driving term with A2 ^ globaUy causes with 
AT[t) ^ a dynamical symmetry breaking. At low, adi- 
abatic driving frequencies w — > 0, we obtain a finite heat 
flux J ^ which is solely induced by dynamical symme- 
try breaking caused by harmonic mixing. Again, the heat 
flux J vanishes in the high frequency limit as expected. 



Adiabatic estimate 

An adiabatic analysis within a linear transport mech- 
anism is no longer possible here: By noting that the 
average temperature -|- Tii{t)]/2 = Tq is time- 

independent, the previous adiabatic estimate now reduces 
to: 



Ja 



to 

UJ 

0; 



2 V- fTL{t)+TR{t)\TL{t)-TRit) 



2n/u 



dt K (To) 



2 

AT{t) 
N 



N 



(13) 



i.e., it predicts a vanishing rachet heat flux. As men- 
tioned above, the observed finite heat flux is due to a 
nonlinear interplay between the driving frequencies for 
which we are even not able to predict a priori the sign. 



Local temperature profile 

To gain further insight into this intriguing regime of 
finite ratchet heat flux, we investigate the local tempera- 
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FIG. 5: (Color-online) Effective temperature profiles of three 
selected frequencies for A2 — 0.3 in Fig. |4]with To = 0.5. 



ture variations across the chain at three distinct different 
driving frequencies tu. The local (kinetic) temperature is 
defined via the equipartition theorem as the time average 
of kinetic energy Teffii) = qf, as in Ref. The effec- 
tive temperature profiles of three numerical runs, denoted 
as (a, b, c) in Fig. IH are depicted in Fig. [5] versus the 
relative site position i/N. In clear contrast to the non- 
driven case with no temperature modulation, a distinct 
temperature profile emerges for the driven case. This av- 
eraged effective temperature T^ffli) lies typically above 
the time-independent average temperature Tav = To with 
regimes of both, positive- and negative-valued local gra- 
dient. Even for the case of slow, adiabatic driving, i.e., 
case (a) in Fig. [51 we cannot now even detect a clear-cut 
mechanism to yield the now positive sign for the shuttled 
heat flux. This corroborates with our reasoning that no 
simple adiabatic estimate can be devised in this situation. 
We also note that with identical signs for Ai , A2 of 



TL(R){t) in Eq. implying that AT{t) = identically, 
no ratchet heat flux can be detected (not shown). 



B. Shuttling heat across a Fermi-Pasta-UIam afi 
lattice 



We next consider the Fermi-Pasta-Ulam a/3 (FPU-a/3) 
lattice, where the potential term V{xi^i, Xi) assumes the 
following form: 



V{Xi+i , Xi) = - (Xi+i -Xi- of 



+ ^(.T,;+i - Xi - aY + ^{xt+i - Xi - a)'^ (14) 



In the simulations, we will employ the dimensionless pa- 
rameters fc = l,Q'=l,/3=l/4. 

The net heat flux J for this system is depicted in Fig. 
ini We again cannot detect any finite heat flux for ^2 = 0. 
For A2 = 0.3, we now detect a negative ratchet heat flux 
in the low frequency, adiabatic limit. This ratchet heat 
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FIG. 6: (Color-online) Ratcheting heat flux in a FPU-a/9 lat- 
tice. Cycle averaged heat flux J as the function of the driving 
frequency uj. The solid circles are the results for symmet- 
ric harmonic driving with A2 — 0. The hollow circles are 
the results for harmonic mixing with A2 = 0.3 and a relative 
phase shift ip = 0. The harmonic driving parameter is set for 
both cases at Ai = 0.6. The simulations are performed on 
a FPU-a/3 lattice with A'' = 50 atoms and with a reference 
temperature To = 1, see text. The arrows indicate the driving 
frequencies used in Fig. [T] 
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FIG. 7: (Color-online) Effective temperature profiles for the 
four selected frequencies with A2 — 0.3 depicted in Fig. [S] 
Here the average bath temperature is held at To = 1. 



flux vanishes as expected in the high frequency limit. 
Surprisingly however, the heat flux does not approach 
zero monotonically as in the case with the FPU-/? lattice. 
At sonic intermediate frequency range, the direction of 
the net heat flux reverses sign! We detect two distinct 
peaks for the heat flux of opposite directions that occur 
within a narrow frequency window. 

We as well depict a plot for the local temperature vari- 
ations across the chain at different driving frequencies uj. 
The effective temperature profiles of four numerical runs, 
denoted as (a, 6, c, d) in Fig. [6l are shown in Fig. [7] ver- 
sus the relative site positions i/N. For frequencies near 
the flux-reversal point, the inner structure of effective 
temperature becomes very complicated and we could not 








FIG. 8: (Color-online) Size dependence of directed heat flux. 
Ratchet heat flux J vs. the driving frequency ui for different 
FPU-a/? lattice sizes iV = 50 (circle), = 100 (triangle) 
and = 200 (square). The driving amplitudes are set at 
Ai = 0.6, A2 — 0.3 and — 0. The estimate of the frequency 
scale set by the thermal response time which scales with the 
anomalous FPU-heat conductivity n{N). 



detect a clear-cut connection between local temperature 
and the direction of ratcheting heat flux. 

When the lattice size is increased, the thermal 
bias AT{t)/N is reduced. One therefore would expect 
a smaller heat flux J. In Fig. [SI we plot J versus lu for 
three different lattice sizes A^ = 50, 100, 200. Contrary 
to common intuition, however, the flnite heat flux J in 
the adiabatic limit lj — > is practically independent of 
system size. This feature again corroborates with the 
fact that in this regime no obvious adiabatic estimate is 
deducible. 

This result of a size-inscnsitivity will prove advanta- 
geous for the experimenters attempting to measure such 
directed ratchet heat flux, e.g. by use of a nanotube cou- 
pled in between two heat contacts. 

The peaks of opposite directions are also found for 
larger system sizes. The positions of these peaks be- 
come red-shifted for larger system size A''. This effect is 
related to the thermal response time as we have detailed 
before with Ref. [l3|. The characteristic frequency can 
be estimated as coc ~ 8ttk{N)/cN'^. The specific heat 
can be approximated as c w 1. The FPU lattice is well 
known to exhibit anomalous heat conduction with size- 
dependent thermal conductivities k{N) [13, [H,!!^]. The 
thermal conductivity k(A^) at temperature T = 1 has the 
following numerical values, k(50) = 19.1,k(100) = 25.1 
and k(200) = 32.9. Thus, the characteristic frequencies 
can be estimated as ujc ~ 0.19 for A^ = 50, ujc ~ 0.063 for 
A^ = 100 and uJc ~ 0.021 for A^ = 200. These values are 
marked as arrows pointing to x-axis in Fig. [51 these val- 
ues impressively corroborate with numerical simulation 
results. 
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FIG. 9: (Color-online) Ratcheting heat flux across a Lennard- 
Jones lattice. Cycle averaged heat flux J vs. the harmonic 
mixing driving frequency u} for a LJ-lattice. Circles are the 
results for the lattice size A*' = 50 in absence of second har- 
monic drive, i.e. A2 = 0. Triangles and squares are the results 
for A'^ = 50 and A*' — 100 with the second harmonic driving 
set at A2 — 0.3. The fundamental driving amplitude is set at 
Ai = 0.6 for all cases. The reference temperature is chosen as 
To = 3 which corresponds to ~ 360 [K]. The arrow indicates 
the driving frequency used in Fig. [8] 



C. Directing lieat across a Lennard-Jones lattice 



We finally consider the physically realistic case of a 
Lennard-Jones (LJ) lattice. The interaction V{xi^i,Xi) 
for a LJ lattice interaction takes the form 



Xi+1 



(15) 

where a is the lattice constant and e is the depth of the 
potential well. New dimensionless parameters can be in- 
troduced by measuring positions in units of a, momenta 
in units of [(me)^/^], temperature in units of [e/ks], 
spring constants in units of [e/a^]. frequencies in units 
of [{e/ma^y^^] and energies in units of e. A particular 
material is described by the pair of parameters a and e. 

The resulting ratchet heat flux J for LJ-lattice is de- 
picted in Fig. [HI Just as is the case in FPU lattice, we 
cannot detect finite flux J for pure harmonic driving with 
A2 ~ 0. In the adiabatic limit a; — > 0, a virtually size- 
independent finite heat flux J results when A2 ^ 0; i.e. 
the emerging ratchet heat flux is rather robust. In our 
simulations we used a dimensionless reference tempera- 
ture To = 3. To give a example, for Argon atoms with 
parameters (a ~ 3A[A], e = 119.8fcB[/\]), this dimension- 
less temperature Tq corresponds to a physical tempera- 
ture T = Toe/kB = 359.4[i4:]. 



1^ 



1^ 
o 



1^ 




(p [units of 27r] 



FIG. 10: (Color-online) Ratcheting heat flux across a 
Lennard-Jones lattice, (a) Cycle averaged heat flux J vs. 
the fundamental driving amplitude Ai with flxed A2 — 0.3 
and ip — 0. (b) Cycle averaged heat flux J vs. the driving 
amplitude A2 with flxed Ai — 0.6 and ifi = 0. (c) Cycle av- 
eraged heat flux J vs. the phase shift Lp with flxed Ai — 0.6 
and A2 = 0.3. For all cases, the lattice size is set as = 50 
and the reference temperature is chosen as To — 3. 



Controlling tieat flux via tfie ptiase stiifl cp 

For single particle Brownian motors, it is well known 
that the directed transport in leading order of nonlinear- 
ity is proportional to the non-vanishing time-averaged 
third moment [1^, H^. Whether this result still holds 
true for a spatially extended ratchet system can be tested 
in the present context numerically only. We conjecture 
that the directed current still will be proportional to the 
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third order moment {/ST {t) / 2TqY , i.e. J oc A\A2COS(p. 
The numerically evaluated heat flux J is depicted vs. 
Ai, A2 and ip in Fig. [101 As expected from theory, it is 
only for small values of driving amplitudes Ai {Ai < 0.5 
in panel (a)) and A2 {A2 < 0.3 in panel (b)) that the 
heat current response J follows the theoretical scaling 
law J oc AfA2 with good accuracy. Note that the sign of 
the flux remains, however, undetermined from the form 
of harmonic mixing with opposite signs for ±A2. At 
larger driving strengths, higher order nonlinear contribu- 
tions yield a sizable contribution, thus causing deviations 
from the leading scaling behavior. Interestingly enough, 
however, the heat current J is found to exhibit the de- 
pendence on the relative phase shift J oc cos if very ac- 
curately, even for substantial large driving strengths; i.e. 
for Ai = 0.6 see in Fig. [TUT c). This numerical finding is 
advantageous for the control of heat current in spatially 
extended systems: the direction of directed heat flow can 
be reversed by merely adjusting the relative phase shift 
if of the second harmonic drive. 



V. 



CONCLUSIONS 



In conclusion, we have demonstrated sizable shuttling 
of a net heat flux across ID homogenous nonlinear lat- 
tices of the FPU-type and Lennard- Jones type by apply- 
ing two different symmetry breaking mechanisms. These 
different mechanisms are imposed via temporal modu- 
lations of the bath temperature(s). For a modulation of 
temperature of one heat bath only a symmetric harmonic 
driving is sufficient to induce a non-vanishing heat flux. 
The resulting ratchet heat flux can be elucidated at low 
driving frequencies by virtue of an adiabatic estimate in 
terms of a single quadrature, see Eq. The expression 
involves the knowledge of the nonlinear temperature de- 
pendent heat conductivity. According to this adiabatic 
analysis, in the linear transport regime, the ratchet heat 
flux is found to be proportional to the square of driving 



amplitude and does vanish in the thermodynamic limit 
iV -> 00. 

For a situation with temperature modulations applied 
at both heat baths, a simple harmonic driving no longer 
induces a ratchet heat flux. If a second harmonic driv- 
ing is included one can break the symmetry dynamically. 
The resulting ratchet heat flux obeys, however, no sim- 
ple adiabatic estimate so that even the direction of the 
ratcheting heat flux cannot be predicted a priori. More- 
over, we detect in the non-adiabatic regime a distinct 
reversal of heat flux for the case of the FPU-a/? chain: 
It occurs around the thermal response time set by the 
anomalous heat conductivity. For the realistic situation 
with a Lennard-Jones lattice we find that the resulting 
flux of the Brownian heat motor is robust and is practi- 
cally independent of system size. 

Noteworthy is the fact that the directed heat fiux is 
substantially larger for a physically realistic Lennard- 
Jones lattice as compared to the situation of coupled, 
non-identical Frenkel-Kontorova lattices, see in Ref. 171. 



As a consequence, an experimental setup as put forward 
with this work seems more feasible to realize such a 
Brownian motor for shuttling heat as compared to a 
physical situation with two coupled Frenkel-Kontorova 
lattices. 
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